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Abstract A sample of 55 distant {z > 0.4) luminous infrared galaxies (LIRGs) selected 
from ISOCAM deep survey fields (CFRS, UDSR, UDSF) have been studied 
on the basis of their high-quality optical spectra from VLT/F0RS2. We em- 
phasize that the only way to derive 0/H metal abundances with sufficient accu- 
racy is by deriving the extinction and underlying B aimer absorption properly, 
on the basis of good-S/N spectra of moderate resolution. Here, the extinction 
coefficient is estimated via two independent methods, e.g., the Balmer-line ratio 
[v4v(Balmer)] and the energy balance between the infrared and H/3 luminosities 
[Av(IR)], are internally consistent, with a median value of 2.36. Oxygen abun- 
dances [12-l-log(0/H)] in the interstellar medium of the sample galaxies have 
been estimated from the extinction corrected emission-line ratios, and show a 
range from 8.36 to 8.93, with a median value of 8.67, which is 0.5 lower than 
that of the local bright disks (i.e., L*) at the given magnitude. A significant 
fraction of distant large disks are indeed LIRGs. Such massive disks could have 
formed ~ 50% of their metals and stellar masses since 2 ~ 1. 

1. Introduction 

The cosmic infrared background resolved by ISOCAM shows that the co- 
moving density of infrared light due to the luminous infrared galaxies (Ljr > 
10^^ Lq; LIRGs) was more than 40 times greater at z ~ 1 than today (Elbaz 
et al. 2002). The main driver for this evolution is the emergence of luminous 
infrared starburst galaxies seen by ISO at z > 0.4 (Flores et al. 1999). Source 
counts and star-formation rates (SFRs) have been studied for the ISO-detected 
objects in the z > 0.4 Universe (Elbaz et al. 1999, Aussel et al. 1999, Flores et 
al. 1999, Franceschini et al. 2003). However, very little was known about the 
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chemical properties of distant LIRGs. In this paper, we present the chemical 
abundances of a large sample of distant (z > 0.4) LIRGs selected from ISO- 
CAM deep survey fields (CFRS, UDSR, UDSF). We estimate their oxygen 
abundances in the interstellar medium (ISM) on the basis of their high-quality 
optical spectra from VLT/FORS2. In particular, this is the first time that the 
luminosity-metallicity (L-Z) relations of such a large sample of distant LIRGs 
are obtained. 

In the local Universe, metallicity correlates well with the absolute luminos- 
ity (stellar mass) of galaxies over a wide magnitude range (e.g., 7 — 9 mag; 
Zaritsky et al. 1994, Richer & McCall 1995, Telles & Terlevich 1997, Contini 
et al. 2002, Melbourne & Salzer 2002, Lamareille et al. 2004, Tremonti et 
al. 2004). In the intermediate-redshift Universe, L-Z relations have also been 
obtained for some sample galaxies. Kobulnicky & Zaritsky (1999) found that 
the L-Z relation of 14 emission-line galaxies with 0.1 < z < 0.5 is consistent 
with that of the local spiral and irregular galaxies. The 16 CFRS galaxies at 
z ~ 0.2 studied by Liang et al. (2004a) fall well into the region occupied by 
the local spiral galaxies. Kobulnicky et al. (2003) obtained the L-Z relation of 
64 intermediate-z galaxies from the Deep Groth Strip Survey (DGSS). In the 
0.6 — 0.82 redshift bin, their galaxies are brighter by ^ 2.4 mag, compared 
to the local (z < 0.1) field galaxies (Kennicutt 1992a,b, hereafter K92a,b; 
Jansen et al. 2000a,b, hereafter J20a,b). This result was confirmed by Maier 
et al. (2004). Kobulnicky & Kewley (2004) obtained the L-Z relation for 204 
emission-line galaxies in the GOODS-N field, and showed a decrease in the 
average oxygen abundance of 0.14 dex from z = to z = 1 for the galax- 
ies with Mb from —18.5 to —21.5. These studies contrast with the result of 
Lilly et al. (2003), who found that the L-Z relation of most of their 66 CFRS 
galaxies with 0.5 < z < 1 is similar to that of the local galaxies from J20b. 
However, Lilly et al. (2003) assumed a constant Ay = 1 in order to account 
for dust extinction. In this study, we investigate the L-Z relation for LIRGs 
at z > 0.4 deteced by ISO, after taking into account their underlying stellar 
absorption and dust extinction properties in detail. 

Throughout this paper, a cosmological model with Hq = 70 km s~^ Mpc~^, 
= 0.3 and il\ = 0.7 has been adopted. Mb is given in the AB magnitude 
system. 

2. Sample selection, observations and data reduction 

The sample galaxies were selected from three ISO deep survey fields: the 
CFRS 3h field, the Ultra Deep Survey Rosat (UDSR) and the Ultra Deep Sur- 
vey FIRBACK (UDSF) fields. In the CFRS 3h field, 70 sources were detected, 
with their 15 /im fluxes in the range of 170 — 2100 /t/Jy (Flores et al. 2004, 
2005, in prep.). The UDSR field refers to the Marano field, which is a deep 
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ROSAT field. FIRBACK is a deep survey conducted with the ISOPHOT in- 
strument aboard ISO, at an effective wavelength of 175 jim. For the UDSR and 
UDSF fields, very deep ISOCAM follow-up observations have been obtained 
(Elbaz et al. 2005, in prep.) reaching flux limits three times lower than for the 
CFRS field. 

In total, 105 objects were selected for our VLT/F0RS2 spectral observations 
with R600 and 1600, at a resolution of 5 A and covering 5000 - 9200 A. The 
slit width was 1.2 arcsec. Spectra were extracted and wavelength calibrated 
using the IRAF package. Flux caUbration was done using 15-minute exposures 
of 3 photometric standard stars per field. 

The redshift distribution of the sample galaxies shows a median value of 
0.587, which is consistent with the results in some other ISOCAM survey fields 
(Flores et al. 1999, Aussel et al. 1999, Franceschini et al. 2003). Fifty five of 
the redshift-identified objects are distant LIRGs with z > 0.4. 

The IR luminosities (and deduced SFRs) were calculated using the proce- 
dure given in Elbaz et al. (2002). The inferred IR luminosity (8 — 1000 fim) of 
the ISOCAM 15/im-detected objects with z > 0.4 shows a similar distribution 
as the local IRAS sample from Veilleux et al. (1995) and Kim et al. (1995), 
with a median value of log (L1R/L0) = 11.32. 

3. A robust estimate of the extinction coefficient 

The extinction inside a galaxy can be derived using its B aimer decrement. 
For the 2; > 0.4 galaxies. Ha is shifted to the near infrared, so that H7/H/3 
can be used to estimate the dust extinction. The stellar absorption under the 
Balmer lines was estimated from synthesized stellar spectra obtained using the 
stellar spectral library of Jacoby et al. (1984). We adopt the interstellar extinc- 
tion law of Fitzpatrick (1999) with R = 3.1, and Case B recombination, with 
a density of 100 cm~^ and a temperature of 10,000 K to estimate the dust ex- 
tinction. Subsequently, the extinction-corrected Balmer emission lines (either 
H/3 or Ha) were used to estimate the SFRs by adopting the calibrations from 
Kennicutt (1998) based on a Salpeter (1955) initial mass function (IMF; O.I 
and 100 Mq mass cut-offs). The derived SFRBaimcr could be compared with 
the SFRiR obtained from the infrared flux. The two SFRs are consistent, and 
their median values are 28 and 31 Mq yr~^, respectively. This kind of consis- 
tency was also confirmed by Hopkins et al. (2003) for the SDSS galaxies, and 
by Kewley et al. (2002) for the Nearby Field Galaxies Survey sample of J20b. 

Because of the large uncertainties related to the measurements of the H7 
line, we need to verify the quality of our derived extinction values. This can 
be done assuming that the infrared data provide a robust SFR estimate for IR- 
luminous galaxies (Elbaz et al. 2002, Flores et al. 2004). We estimate a new 
dust extinction coefficient, ylv(IR), by comparing SFRjr and SFR2.87H/35 the 
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energy balance between the IR and H/3 luminosities. Figure 1 shows that the 
derived Ay (IR) is consistent with Av(Balmer) for most galaxies, most of them 
falling in the ±0.64 rms discrepancy. The derived median value of Av'(IR) is 
2.36 for the z > 0.4 galaxies. 




Figure 1. The relation between the extinction values derived from the Balmer decrement 
[^v(B aimer)], and from the energy balance between the IR radiation and the optical H/3 
emission-line luminosities [Av(IR)]. They are consistent with each other. The two dashed lines 
refer to the results with ±0.64 rms. 



4. Abundances in the interstellar medium and the 
luminosity-metallicity relation 

The chemical properties of the gas and stars in a galaxy are like a fossil 
record chronicling its history of star formation and its present evolutionary 
status. The high-quality optical spectra from VLT/F0RS2 make it possible to 
obtain the chemical abundances of the ISM for these distant LIRGs. This will 
be the first data set of chemical abundances of a large sample of distant LIRGs. 

The oxygen abundances of the distant LIRGs were estimated via the "strong- 
line" method (i?23 and O32), using the calibration of Kolbinicky et al. (1999) 
for the metal-rich branch. The derived 12-i-log(0/H) values of the galaxies 
range from 8.36 to 8.93, with a median value of 8.67. 

Figure 2a compares the L-Z relation for the LIRGs to that of local disks 
(K92b, J20b), which are restricted to moderately star-forming [EW(H/3) < 20 
A] galaxies, following Kobulnicky et al. (2003). Here we did not include 
other local samples selected based on their UV or Ha emission; they mostly 
include low-luminosity (low-mass?) systems in the local Universe (Contini et 
al. 2002, Lamareille et al. 2004, Melbourne & Salzer 2002). We could not 
compare directly with the large sample of SDSS galaxies given by Tremonti et 
al. (2004) because of the large scatter of the data. 
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The distant LIRGs exhibit ~ 0.3 dex (~ 50%) lower oxygen abundances 
than local star-forming galaxies at the given magnitude (the median value). 
P^GASe2 models (Fioc & Rocca-Volmerange 1999) predict a total mass rang- 
ing from 10^^ Mq to < 10^^ M0 for the LIRGs, which can be twice the stellar 
masses of distant LIRGs (1.4 x 10^° - 2.9 x 10^^ M0) derived by Zheng et 
al. (2004) on the basis of i^-band luminosities. About 36% of LIRGs are large 
disks; a similar fraction (about 32%) is estimated from the sample of large 
disks of Lilly et al. (1998). These massive LIRGs have high SFRs. The time- 
scale to double their stellar masses Tsf can be short, i.e., 0.1 — 1 Gyr. Such 
massive disks could have formed ~ 50% of their metals and stellar masses 
since 2; ~ 1. Hammer et al. (2004, 2005) have investigated whether the LIRG 
properties could be related to recent and significant star formation in massive 
galaxies, including spirals. 

These distant LIRGs were also compared with two samples of galaxies in a 
similar redshift range, taken from Kobulnicky et al. (2003) [EW(H/3) < 20 A] 
and Lilly et al. (2003). Kobulnicky et al. (2003) estimated the O/H values 
using the i?23 and O32 parameters obtained from the corresponding equivalent 
widths of the emission lines, which are beheved to be less affected by dust 
extinction (Kobulnicky & Phillips 2003). The metallicities of their galaxies 
are similar to ours, but the galaxies are fainter at a given metallicity. This 
discrepancy in Mb reflects that our sample galaxies are brighter and possibly 
more massive than the rest-frame blue-selected sample of DGSS galaxies (see 
Liang et al. 2004b for details). 

The comparison with Lilly's sample helps us to understand the strong effects 
of dust extinction on the derived oxygen abundances of galaxies. In Fig. 2b, 
Lilly's sample has been restricted the galaxies in the CFRS 3h and I4h fields, 
the two fields surveyed by ISOCAM. The 10 ISO galaxies (the solid trian- 
gles) among the 42 sample galaxies show a median value of 12-i-log(0/H) = 
8.98, which is ~ 0.3 dex higher than the median value of our distant LIRGs. 
Indeed, Lilly et al. (2003) assumed a constant extinction of Ay = 1 for 
all their galaxies, which can underestimate the average extinction for LIRGs 
(which has a median value of 2.36), and consequently leads to underestimated 
[Oii]A3727/H/? ratios, and hence to overestimated oxygen abundances. To 
estimate the chemical abundances of such distant LIRGs, underlying stellar 
absorption of Balmer lines and dust extinction should be considered carefully. 
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Figure 2. The Ms-metallicity relation of our distant LIRGs (with a typical uncertainty of 
0.08 dex on metallicity), compared with other samples: (a) with the local galaxies from K92b 
and J20b; PJ;gase2 infall models are superimposed, assuming a total mass of 10^^ Mq and 
infall times of 5 Gyr and 1 Gyr (solid and dashed lines with pentagons, respectively); (b) with 
the galaxies of Lilly et al. (2003). The linear least-squares fits to the samples are also given. 
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